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Abstract 

We extend our study of bottomonium spectral functions in the quark- 
gluon plasma to nonzero momentum. We use lattice QCD simulations 
with two flavours of light quark on highly anisotropic lattices and treat the 
bottom quark with nonrelativistic QCD (NRQCD). We focus on S wave 
(T and channels and consider nonrelativistic velocities, v/c < 0.2. A 
comparison with predictions from effective field theory is made. 
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1 Introduction 



Quarkonia, heavy quark-ant iquark bound states, play an important role as probes 
of the quark-gluon plasma (QGP) created in relativistic heavy-ion collisions at 
RHIC and the LHC pQ (for reviews, see e.g. refs. [2j[3]). How the width of 
a particular bound state - and its dissociation rate - varies with temperature, 
depends on the quark content (charmonium, bottomomium) and on its quantum 
numbers (e.g. S or P wave, in nonrelativistic notation). Recently the CMS 
collaboration has observed sequential Upsilon suppression in PbPb collisions at 
the LHC [HE], which has inspired significant phenomenological and theoretical 
activity [6HTTj. 

Traditionally, quarkonium suppression has been studied using potential mod- 
els (see refs. [12], [13] and references therein) and with lattice QCD computa- 
tions of quarkonium spectral functions [HH2Q] (recent reviews can be found in 
refs. [21] [22]). Several years ago the theoretical study of quarkonia in a ther- 
mal medium was formulated in a more systematic fashion by casting the prob- 
lem in the language of effective field theory (EFT) [23H31] . Besides the energy 
scales available in vacuum: the heavy quark mass M q , the inverse system size 
l/a ~ M q a s , and the binding energy E b ~ M q a 2 s (with a s the strong coupling 
constant), new thermal scales are provided by the temperature T and the inverse 
Debye screening length 1/rp ~ ^/a s T. The relevance of these scales depends on 
the temperature of the QGP and the magnitude of the strong coupling constant. 
In analytical studies, a s is usually assumed to be small enough for a hierarchy 
of scales and EFTs to emerge. An important conceptual finding was the emer- 
gence of a complex potential at nonzero temperature [23], which has stimulated 
further complex potential model studies [32j[33] as well as attempts to extract 
the complex potential from lattice QCD [M[l35]. 

For temperatures reached in heavy-ion collisions and in the case of bottomo- 
nium, the heavy quark mass is the highest energy scale present. Integrating out 
this scale yields nonrelativistic QCD (NRQCD) [36ti39j . just as in vacuum. Re- 
cently we have employed lattice NRQCD [^014^2] to study the fate of P wave [13] 
and S wave [I3HH] bottomonium states nonperturbatively, using lattice QCD 
simulations of nonrelativistic bottom quarks, propagating through a medium of 
thermal gluons and two light flavours, at temperatures between 0.4T C and 2.1T C . 
We found that the use of NRQCD greatly enhances the signature for quarkonium 
melting/survival, since it avoids several problems which have complicated the 
study of relativistic quarks in thermal equilibrium [4"5l - f4"8] . 

Our studies indicate that bound states in the P wave channels (xb, hf,) melt 
quickly as the temperature is raised above T c , while the signal in the S wave chan- 
nels (T,?]^) is consistent with survival [43J. A closer study of the S wave spectral 
functions, constructed using the maximum entropy method (MEM) [IH1I50], sub- 
sequently showed that the ground state peaks appear to survive, whereas excited 
states are suppressed, consistent with the recent CMS results (HE]. Moreover, 
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medium effects in the ground state peak can be captured by a temperature- 
dependent position and width jH]. The thermal mass shift and width were 
found to be consistent with those predicted by EFT calculations [27], assuming 
a s ~ 0.4 |H]. 

In this paper we extend the lattice calculation and consider bottomonium 5* 
wave correlators at nonzero momentum. While we find a significant momentum 
dependence in the euclidean correlators (as expected), we observe that this mo- 
mentum dependence is to a large extent independent of the temperature. This 
result is further investigated by constructing the corresponding spectral func- 
tions at nonzero momentum and extracting the momentum-dependent position 
and width of the ground state peak. We note here that in the literature only 
a handful of lattice studies of spectral functions at nonzero momentum can be 
found, all using the relativistic formulation: see ref. [51] for light (staggered) 
quarks and refs. [52JE3] for studies of charmonium. A lattice study of the (gauge 
fixed) quark propagator at nonzero momentum can be found in ref. [5*4] . 

This paper is organised as follows. In the next section, we first discuss how 
the different scales appear in our lattice simulation and, most importantly, what 
ground state velocities can be reached. The main results for the correlators and 
spectral functions are given in section [3] for the hadronic phase and in section H] 
for the quark-gluon plasma phase, where we also comment on possible compar- 
isons with the predictions from thermal EFT calculations at nonzero velocity. A 
summary is given in section |5j 

2 Scales on the lattice 

In this section, some simple estimates of the scales appearing in the lattice sim- 
ulation and in the EFT formulation are given. We consider bottomonium: for 
the sake of simplicity let us assume a quark mass M q ~ 5 GeV and an S wave 
groundstate (at rest) of Ms ~ 9.5 GeV. The temperatures that we are able to 
reach in our current two-flavour, highly anisotropic lattice simulations go up to 
458 MeV. Some details of the lattice setup are given in table [TJ The temperature 
is very precisely determined in MeV; in units of T c less so, primarily because we 
only have a rough estimate of T c on these lattices. More details can be found in 
refs. [T8" | l44 p52ll5"5] . At these temperatures, one may use 0.3 < a s {T) < 0.4 (previ- 
ously we found a s ~ 0.4 [H], based on a comparison with EFT predictions [27]). 
We arrive therefore at the hierarchy 



where the final comparison depends on the size of a s and on the temperature. The 
use of NRQCD is therefore very well motivated. Details of the 0{v A ) improved 
lattice NRQCD formulation we use can be found in ref. |4"4"] . 



M q > M q a s > T ~ M q a 2 s , 
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Table 1: Two-flavour lattice details: the lattice size is x iV T with N s = 12, 
and the lattice spacing is a s ~ 0.162 fm, a™ 1 = 7.35(3) GeV determined from 
the IP — IS* splitting in charmonium. The anisotropy is a s /a T = 6 [H]. The 
temperatures in MeV have uncertainties of ~ 0.5% from the uncertainty in the 
lattice spacing, while the temperatures in units of T c are rough estimates with 
an accuracy of ~ 10%. 
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Table 2: Nonzero momenta used in this study. Also indicated are the corre- 
sponding velocities v = \p\/M s of the ground states in the vector (T) and pseu- 
doscalar (rj b ) channels, using the ground state masses determined previously [43] . 
My = 9.460 GeV and ikL, = 9.438 GeV. 



The momenta and velocities that are accessible on the lattice are constrained 
by the discretization and the spatial lattice spacing. The lattice dispersion rela- 
tion reads 

22 , • 2 Pi 27m, N s N s 

cCp = 4 > sin — , pi = — , <rii< — . 2.2) 

s r £^ 2 , Ft N ■ , 2 _ 2 v ) 

i=l 

To avoid lattice artefacts, only momenta with rij < N s /4: are used: we consider the 
combinations (and permutations thereof) given in table EJ The largest momen- 
tum, using n = (2,2,0), is |p| ~ 1.73 GeV, corresponding to v = \p\/Ms — 0.2. 
Therefore, the range of velocities we consider is non-relativistic. 

3 Low temperature 

We start by presenting the results at the lowest temperature in the hadronic 
phase, i.e. T/T c ~ 0.42. Fig. [T] (left) shows the euclidean correlation functions 
at nonzero momentum, normalized by the correlator at zero momentum. The 
momentum dependence is introduced in such a way that G(t = 0, p) = G{r = 
0,0) for all momenta. 
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Figure 1: Euclidean correlators at nonzero momentum p = |p| normalized by 
the correlator at zero momentum, G(t, p)/G(r, 0), at the lowest temperature in 
the vector (T) channel (left) and the pseudoscalar (77;,) channel (right). Error 
bars are considerably smaller than the symbols. 



We observe considerable momentum dependence. This is easily understood. 
At large enough euclidean times, the correlator is dominated by the ground state 
and 

G(r,p) ~ e - M ^ T . (3.1) 



We then find 



where 



G(t, p) 
G(t,0) 



e 



-AM(p)r ^ x _ AM ( p ) r + . . . j (3.2) 



AM(p) ee M(p) - M(0) » = V^ 2 . (3.3) 

Here Mki n is the kinetic mass, which is equal to the rest mass M(0) when the 
quark mass is carefully tunedjj The ratios shown in fig. [1] are therefore expected 
to drop to zero exponentially at large euclidean time. However, at a finite time 
and for small momenta, i.e. when AM(p)r <C 1, this translates into linear time 
dependence, as indicated in eq. (13. 2p and visible in fig. [T] as well. 

In NRQCD, euclidean correlation functions and spectral functions are related 

by 

/•"max J,, 

G(t, p) = / — K(t, u>)p(u>, p), K(t, co) = , (3.4) 

both at zero and nonzero temperature [25 l l4~3 l l4"4"]. We construct the spectral 
functions using the maximum entropy method jlHl[50], referring to ref. [S] for 
details of MEM in this context. We emphasize that in this setup the temperature 
dependence does not enter via the kernel K(t, oj), but arises solely from the pres- 
ence of the medium of gluons and light quarks at different temperatures. This 



1 In our study this holds only approximately. 
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Figure 2: Spectral functions p(oj,p), normalized by the heavy quark mass, 
as a function of energy in the vector (T) channel at the lowest temperature, 
for several momenta p. The graph on the right shows a close-up, with vertical 
lines indicating the position of the ground state at each momentum obtained via 
standard exponential fits. 
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Figure 3: As in the preceding figure, for the pseudoscalar (rjb) channel. 

greatly enhances the robustness of MEM [H] and avoids a number of problems 
associated with relativistic quarks, such as zero- modes [16]. Note that due to 
the separation of scales (M q 3> T), thermal effects for the heavy quark are expo- 
nentially suppressed; such terms are in fact removed from the NRQCD partition 
function when integrating out the heavy quark energy scale. Whether b quarks 
actually thermalise in the quark-gluon plasma created in heavy-ion collisions at 
RHIC and LHC is a separate question which we will not address here. 

The spectral functions are shown in figs. [2] and [3] for the vector and the 
pseudoscalar channel. We observe that the ground state peak shifts to the right 
with increasing momentum, as expected. The vertical dashed lines in the figures 
on the right indicate the position of the ground state at each momentum obtained 
via standard exponential fits. Agreement between both methods can be seen. 
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The second bump in the figures on the left can be identified with the first excited 
state [H]. A detailed analysis of systematic uncertainties associated with MEM 
can be found in ref. |44j . 



4 High temperature 



We now turn to temperatures in the quark-gluon plasma phase. Here we find that 
the correlators depend on both temperature and momentum. This is illustrated 
for the vector channel in fig. HJ on the left we present the temperature depen- 
dence via the ratio G(t, p; T)/G(r, p; T ) for fixed p = 1.38 GeV and reference 
temperature To = 0.42T C , while on the right we show the momentum dependence 
via the ratio G(r,p;T)/G(r,0;T) for fixed T = 1.68T C . We observe that the 
temperature dependence at nonzero momentum is similar to the one found at 
vanishing momentum (compare with fig. 1 of ref. jH]) and is of the order of a 
few percent. Similarly the momentum dependence in the quark-gluon plasma is 
similar to what we found at low temperature (compare figs. [T]and|l]). 
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Figure 4: High-temperature results in the vector (T) channel. Left: Euclidean 
correlators at high temperature normalized by the correlator at the lowest tem- 
perature (T/T c = 0.42) at fixed momentum p = 1.38 GeV. Right: Euclidean 
correlators at nonzero momentum normalized by the correlator at zero momen- 
tum at fixed temperature T = 1.68T C . 

In order to analyze whether the momentum dependence changes as the tem- 
perature is increased, we consider the following double ratios 

G{r,p;T) / G(T,p;To) 

G(t,0;T)/ G(t,0;T )' 1 " ] 

where the reference temperature is again To = 0.42T C . The results are shown in 
fig. [5] in the vector channel (the pseudoscalar channel is similar) . We observe that 
after the normalization, the remaining momentum dependence is very mild, and 
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always less than 0.5% at the largest euclidean time. This apparent temperature 
independence of the momentum dependence provides a clear prediction which 
can be contrasted with potential models and EFT calculations. 

The corresponding spectral functions are presented in figs. [6] and [71 As ex- 
pected from the results presented above and in ref. [Bj, we find that the ground 
state peak moves to the right as the momentum is increased and that it broadens 
and reduces in height as the temperature is increased. As in ref. jH], the first 
excited state and other features at higher energy quickly become suppressed as 
the system heats up. We note that the ground state peaks are very stable at 
low and intermediate temperatures. This is less so at the highest temperatures, 
presumably due the limited number of temporal lattice points available jH] . 

Following ref. jH], we extract the position M(p) and width T(p) of the ground 
state peak from the spectral functions. The results are shown in fig. [8] in the 
vector and pseudoscalar channels, as a function of v 2 , where v = |p|/M(0,T), 
with M(0, T) the peak position at zero momentum and temperature T. The 
temperature dependence of M(0, T) has been discussed and compared with EFT 
calculations in ref. jH]. The peak position below T c has been extracted using 
standard exponential fits, but agrees with the one obtained from the spectral 
functions, see figs. [2] and Note that the width is normalized by the temperature. 
To avoid cluttering, error bars are indicated for three velocities only; each of 
these has two error bars representing the uncertainty due to the number of time 
slices that are included in the MEM analysis (systematic uncertainty, left error 
bar), and due to the number of configurations used in the MEM computation 
(statistical uncertainty, right error bar); there is no significant dependence on the 
default model, even at the highest temperature. Errors are under control at the 
lower temperatures, but uncertainties become large at the highest temperatures. 
The position of the peak is easier to determine than its width. Ref. [H] contains 
a detailed discussion of error bars and uncertainties. 

We observe that the peak position increases linearly with v 2 , as expected. 
Assuming the lowest-order, nonrelativistic expression M(p) = M(0) + p 2 /2M(0), 
one finds 

M(v) p 2 1 9 . . 

Mf = 1 + WW) = 1+ 2" ' (4 ' 2) 
which is indicated with the dotted lines in the left figures. The slope of the data 
points at the lowest temperature is slightly less than 1/2, which can be improved 
by a careful tuning of the heavy quark mass. Systematic uncertainties are too 
large to find a trend in the velocity dependence of the width; we note that the 
width is independent of the velocity within errors. 

The dependence on the velocity can be compared with EFT predictions. In 
ref. [31] a study of the velocity dependence was carried out in the context of QED, 
working in the rest frame of the bound state (i.e. the heat bath is moving). In 
order to compare with our setup, we consider the case in which the temperature 
is low enough for bound states to be present and that the velocities are nonrela- 
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Figure 5: Normalized euclidean correlators in the vector (T) channel: shown is 
the double ratio [G(r, p; T)/G(r, 0; T)]/[G(r, p; T )/G(t, 0; T )], where the refer- 
ence temperature T = 0.42T C . Each panel shows several momenta p at a fixed 
temperature T. 
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Figure 6: Spectral functions p(u, p), normalized by the heavy quark mass, as a 
function of energy, at the six different temperatures above T c in the vector (T) 
channel, for several momenta. 
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Figure 7: As fig. [61 for the pseudoscalar (r/ b ) channel. 
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Figure 8: Position of the ground state peak M(p)/M(0) (left) and the upper 
limit on the width of the ground state peak, normalized by the temperature, 
T/T (right), as a function of the velocity squared (v 2 ) in the vector (T) channel 
(above) and the pseudoscalar (rj b ) channel (below), as extracted from the spectral 
functions. The ground state positions below T c are obtained using a standard 
exponential fit. The dotted line in the left figure represents M(p)/M(0) = l + 2 . 
Error bars are indicated for three velocities only; each of those has two error bars 
representing statistical and systematic uncertainties jH]. On the left, error bars 
are not shown for the highest temperature; they exceed the graph size. 



tivistic. In that case, one finds [31], in the rest frame of the bound state and at 
leading order in the EFT expansion, 



E^^EZ log (I±^), (4 .3) 

r 2v \i-vj 

where r is the width at rest. Interpreting the width as an inverse lifetime, one 
can express this result in the rest frame of the heat bath by dividing with the 
Lorentz factor 7 = — v 2 . An expansion for nonrelativistic velocities then 

yields 

If we take this result and apply it to our study of bottomonium, we find that the 
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effect of the nonzero velocity shows up as a correction at the percent level (recall 
that v 2 < 0.04), which is beyond our level of precision but consistent with the 
observed v independence within errors. Similarly, additional thermal effects in 
the dispersion relation are currently beyond our level of precision. 

5 Summary and outlook 

In this paper we extended our analysis of bottomonium in the quark-gluon plasma 
to nonzero momentum, using lattice QCD simulations with two flavours of light 
quarks and nonrelativistic dynamics for the bottom quark. We analyzed both 
the euclidean correlators as well as the associated spectral functions constructed 
using the maximum entropy method. While we observed both temperature and 
momentum dependence in the correlators and spectral functions, we found that 
the momentum dependence is effectively temperature independent. This is seen 
directly in the correlators as well as in the position and width of the ground state 
and can be compared with existing and future EFT predictions. 

As an outlook we plan to improve the tuning of the heavy quark mass, with the 
aim to enhance the prospects of detecting thermal deviations from the standard 
dispersion relation. We are planning to carry out this analysis in a quark-gluon 
plasma with Nf = 2 + 1 rather than 2 flavours, with both a smaller spatial lattice 
spacing and a larger spatial extent, which will also allow us to reach higher 
velocities. 

We hope that this work provides further encouragement to study quarko- 
nia at nonzero momentum using EFT and potential model approaches and can 
contribute to the interpretation of the experimental results for bottomonium in 
heavy ion collisions at the LHC. 
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